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Abstract 



D , Various primary user (PU) radios have been allocated into fixed frequency bands in the whole spectrum. 



A cognitive radio network (CRN) should be able to perform the wideband spectrum sensing (WSS) to detect 
temporarily unoccupied frequency bands. We summarize four occupancy features for the frequency bands. 1. The 
occupancy is sparse; 2. The frequency band allocation information is fixed and common; 3. There are three categories 
, for the frequency band usages; 4. The occupied frequency bands are common in the CRN. For the first time, we 

O 

consider all features as the prior knowledge in the compressed sensing based cooperative WSS (CWSS) algorithm 

■ design for a centralized CRN. We propose a modified orthogonal matching pursuit (Mod-OMP) algorithm and a 
I modified simultaneous orthogonal matching pursuit (Mod-SOMP) algorithm for the CWSS. We compare the CWSS 

■ performance of Mod-OMP/Mod-SOMP with the original OMP/SOMP and show the performance improvements. 

(N 
O 

I. Introduction 

According to the US frequency allocation chart [0, a variety of primary user (PU) radios have been 



• allocated to the frequency bands covering the whole spectrum. However, existing PUs are not occupying 
all frequency bands all the time 0. For example, filed measurement of a 2.5 GHz spectrum shows that the 
occupancy rate can be as low as 10% [1 IJ. Cognitive radio network (CRN) is a network of secondary user 
(SU) using the temporalily unoccupied frequency bands. A CRN should be able to perform the wideband 
spectrum sensing (WSS) to detect those unoccupied frequency bands. 

One challenge for the WSS is that the possible unoccupied frequency bands are spread over a multi- 
GHz spectrum. If the WSS is performed following the Nyquist sampling rate, it will be infeasible for 
current analog-to-digital converters (ADC). Several features can be utilized to overcome the sampling 
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rate limitation. Feature 1 is that the spectrum is sparse, which means that the spectrum occupancy rate is 
low. This has made the compressed sensing based algorithms possible for the WSS. Several compressed 
spectrum sensing (CSS) receivers [fTOl . [[T2l have been proposed to perform the WSS under sub-Nyquist 
sampling rate. Feature 2 is that the frequency bands allocated to the PUs are fixed and can be considered 
as the prior knowledge to the SU. Reference BEl has made use of such prior knowledge and proposed 
a basis pursuit (BP) [2] based algorithm, the mixed I2/I1 norm denoising operator algorithm (MNDO), 
reaching a better reconstruction performance than the traditional reconstruction algorithms. However, our 
previous work in [[T6l shows that when a third feature is considered, the reconstruction performance can 
be further improved. Feature 3 is that the usages of the frequency bands can be grouped into 3 categories 
0: 

• Some frequency bands in the spectrum are largely unoccupied most of the time. 

• Some other frequency bands are only partially occupied. 

• The remaining frequency bands are heavily used. 

In [fT6ll , we have demonstrated that feature 3 can be integrated into the classical greedy algorithm, orthog- 
onal matching pursuit (OMP) [fT3l , and we have proposed a modified OMP (Mod-OMP) reconstruction 
algorithm for the single user WSS. 

When multiple SUs can cooperate, the cooperative WSS (CWSS) is a common approach to achieve 
better WSS performance in a CRN. An additional feature needs to be pointed out for the CWSS: Feature 
4, the occupied frequency bands are common among all SUs within the same CRN. The CRN can be 
centralized or decentralized. Here we consider the centralized CRN with a fusion center. All SUs can only 
cooperate with the fusion center to perform the WSS. Moreover, we consider the greedy algorithms due 
to its low computation complexity [fT3l . To the best of our knowledge, the CWSS algorithms using greedy 
algorithms for the centralized CRN are very limited. In [fT4l . a distributed compressed sensing (DCS) 
O based architecture is proposed for the centralized CRN. It uses the simultaneous orthogonal matching 
pursuit (SOMP) reconstruction algorithm to perform the CWSS. However, it only takes advantage of 
feature 1 and feature 4. 

In this correspondence, we focus on the design of the CWSS algorithm using CSS for a centralized CRN. 
For the first time, all four features are integrated. Two common schemes for the CWSS are considered: 
decision fusion and data fusion [15J. For the decision fusion CWSS, we make use of the single user 
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Mod-OMP and propose the Mod-OMP based decision fusion algorithm (Mod-OMP DeF). For the data 
fusion CWSS, we propose a novel modified SOMP (Mod-SOMP) algorithm as well as the Mod-SOMP 
based data fusion algorithm (Mod-SOMP DaF). We adapt the Hanssen-Kuiper skill score or R score flSl 
metric to evaluate the accuracy of the decisions made upon the frequency band occupancy. Simulations are 
conducted to compare the performances of the decision fusion algorithm with original OMP (OMP-DeF), 
the data fusion algorithm with original SOMP (SOMP-DaF), as well as the proposed Mod-OMP DaF and 
the Mod-SOMP DeF. 

This correspondence is organized as follows. In Section Ull we review the system models of the CSS 
and the R score metric. In Section Hill we introduce the four CWSS algorithms in detail. Simulation 
results for all algorithms are shown in Section |IV]and conclusions are made in Section |Vl 

II. System Model 

Fig. [T] depicts an example of a centralized CRN with 1 fusion center and 3 SUs. There are 3 PUs 
activated within the CRN. Each PU can only occupy the frequency band allocated to itself. The SUs in 
the CRN can use any of the unoccupied frequency bands. In the CRN, all SUs and the fusion center 
perform the CWSS to determine the temporarily unoccupied frequency bands. The bottom-right of Fig. \T\ 
shows the spectrum observed at each SU when 3 PUs are activated. Feature 4 can be observed, that the 
occupancy of the frequency bands are common among all the SUs. 

A. Compressed Spectrum Sensing 

One challenge for the WSS is the sampling rate. Recent progress in compressed sensing [|T1, dH provides 
a framework to overcome the sampling rate problem. Here we briefly review the basic CSS models. 

1 ) Single User Compressed Spectrum Sensing Model: Assume there is an analog signal x (t) that can 
be represented by a finite weighted sum of basis functions ^pi (t) with coefficients s (i): 

N 

x{t) = Y.sit)^dt) (1) 
1=1 

In the discrete time framework. Equation [1] can be represented as x = ^'s, where ^ is the N x N matrix 
with each column as the discrete time representation of basis function ipi (t), x is the x 1 vector, the 
discrete time representation of x (t), and s is an x 1 signal vector with A^ entries s (i). After sampling 
X with an M X A^ matrix we have the A^ x 1 measurement vector y: 

y = $X = = 0S (2) 



For convenience, we call © the sampling matrix. If s is sparse, and if meets the restricted isometry 
property (RIP) [[11, H, then s can be reconstructed perfectly with M < N measurements by solving the 
following /i-norm optimization problem: 

s = argmin ||s||^ s.t. y = 0s (3) 

where is the Ip norm defined as: 

ii-ii, = (E(-r)' (4) 

In the CSS, s can be viewed as the spectrum observed at the SU, denotes the sampling method 
and y is the measurements at the receiver. If meets the RIP, s can be reconstructed perfectly via 
Equation [3] with M « N measurements, and therefore sub-Nyquist rate sampling can be achieved. The 
details for the design of can be found in [[TOl , [[T2l . In this correspondence, however, we focus on the 
reconstruction algorithms solving Equation [3] We will integrate the specific prior knowledge from feature 
2 - feature 4 into OMP and introduce the Mod-OMP for decision fusion. 

2) Cooperative Compressed Spectrum Sensing Model: When J signal vectors Sj are sampled with the 
same at each SU, we have the M x J measurement matrix Y and the N x J signal matrix S: 



yi y2 ••• yj 







Si S2 ... Sj 



0S (5) 



Such sampling model is referred to the DCS [[3, [[T4l . Similar to the single user CSS, S can be recon- 
structed from Y with M « N measurements. The joint-sparsity-model 2 (JSM-2) is of special interest, 
where the index set {i\sj (i) ^ 0} are common for all while with different values. Such commonness 
exactly fits feature 4, and SOMP has been adapted to reconstruct S [[T4l . Again, we will focus on the 
reconstruction algorithm and introduce the novel Mod-SOMP for data fusion considering the specific prior 
knowledge provided from feature 2 - feature 4. 

III. Cooperative Wideband Spectrum Sensing Algorithms 

In this section, we describe the representations of the spectrum using feature 2 and summarize feature 
3 in Table [B Then, the R score metric is reviewed. Finally, the Mod-OMP DeF and the Mod-SOMP DaF 
using feature 2 - feature 4 will be introduced respectively. 



A. Representation of the Frequency Allocation Information 

Assume a centralized CRN has J SUs and a fusion center. The B Hz spectrum available for the CRN 
is allocated with K PUs. The discrete spectrum measured at the j*'^ SU can be described using s,: 



Si 



1 T 



(6) 



where T is matrix transpose and sj (i) are the uniformly sampled frequency values over the B Hz spectrum 
According to feature 1, Sj is sparse: 



s 



,„o/^^<10%, j = l,2,---,J (7) 



The frequency allocation information in feature 2 can be represented by index sets u^: 

ui = {l,2,---6i} 
U2 = {61 + 1,61 + 2, ■■■62} 



(8) 



Uk = {bk-i + l,bk-i + 2,---N} 
where 61 ■ ■ ■ 6fc_i are frequency indices determining the boundaries of different frequency bands. 

Feature 3 is summarized in Table |I] via time and frequency, can be grouped into the category sets 

C„ = {ufc \k G Category n} ,n = 1,2,3 

(9) 

^2 = {1,2,-- - ,iV} = UC„ (C,nC, = 0,for^^j) 

n 

In the noiseless case, if \\sj (i)|| > and if i G u^, then frequency band u^. is occupied and all i within 
should not be used by the CRN. Let the K x 1 vector o denote the occupancy decisions of the frequency 
bands. If Uk is occupied, then o (k) = 1, otherwise o (k) = 0. 

If decision fusion is used, each SU will first obtain a local estimate 6^, then the fusion center will 
obtain an overall estimate 6. For data fusion, 6 is obtained at the fusion center directly. 

Since there are K decisions to be made for the whole spectrum, we need to consider the following 
decision statistics [|3l, [|9l from o and 6: 

• Number of Hit, a, the number of correct decisions for o{k) = 6{k) = 1. 

• Number of Miss, 6, the number of wrong decisions for o(A;) = 1, 6(A;) = 0. 

• Number of False, c, the number of wrong decisions for o(A;) = 0, 6(k) = 1. 

• Number of Correct, d, the number of correct decisions for o(A;) = 6(A;) = 0. 



The Hanssen-Kuiper or R score is a common metric using the above statistics to evaluate the perfor- 
mance of correct decisions and wrong decisions. The R score is calculated as: 

R^wi^-W2-^ (10) 
a + c + a 

where wi and W2 are the weights denoting the credit to correct decisions and the penalty to wrong 
decisions, respectively. In this correspondence, we set = W2 = 1, meaning that the loss due to wrong 
decisions is equal to the benefit from correct decisions. Other weights can be set according to desired 
situations. In such setup, it! e [—1,1], with 1 as the score indicating that all decisions are correct, while 
—1 as the score indicating that all decisions are wrong. Next, we will modify OMP and SOMP with the 
prior knowledge and C„. 

B. The Mod-OMP for Decision Fusion 

Mod-OMP is modified upon the original OMP. For the j*'* SU, Mod-OMP starts with local measurement 
vector Yj, sampling matrix 0, error tolerance r] determined by noise level, and the index sets and C^. 
Mod-OMP gets the estimate Sj from the following steps: 

1) Initialize: Set the iteration counter t = 0. Initialize the residual resj*-* = y^, 0^*^ = 0, the index 
set A(*) = Ci. If A(*) = 0, jump to Step 4. 

2) Solve the least-squares problem to obtain z^: 

Zj = argmin 0(*)z - y^- (11) 

Z Z 

3) Update the new residual: 

resf = y^- - Q^^^zj (12) 

4) Increment t. Find the index A^*) that satisfies the following equation: 



A^*^ — arer max 

°i=l,...,N 



resf '\9,) (13) 



where 9i denotes the i*'* column vector of 0. (a, /3) denotes the inner product of two vectors a and 
f3. If the maximum occurs for multiple indices, break the tie deterministically. 

5) If A(*) e Ufe C C2, let = A(*-i) Uufe. If A^ e ujk C C3 or C„ = 0, let A^ = A(*-i) U A^. 

6) Set the matrix of chosen atoms 0^*^ = 0a(*)' where 0a(*) is a sub-matrix of containing columns 
with indices in A^*^ 



7) Obtain Zj using Equation [TT] and update the new residual res^*^ using Equation [121 

8) If t < m or res^*"* ^> return to Step 4. Otherwise, Zj is the estimated signal Sj with non-zeros 
at indices listed in A*^*^ 

Notice that if C„ = and = 0, there is no prior knowledge available, and Mod-OMP becomes 
the original OMR It can be seen that OMP only finds one non-zero index i in each iteration and treats 
all i independently. The modifications made by Mod-OMP using the prior knowledge and C„ are 
described as follows, i G Ci are included in the initial index set as A^^^ = Ci, because values under 
those indices are non-zero all the time. For i E C2, if a single i is selected in iteration t and if i G u^, then 
A(*) = A*^*~^) Uufe, because there must be other non-zeros values under frequency band u^. For i G C3, 
only one i will be selected in each iteration, because frequency bands in C3 are only partially occupied. 

The modifications above have the following advantages: 

1) The number of iterations are reduced from the number of occupied frequency points to the number 
of occupied frequency bands. 

2) The reconstruction errors are restricted within several selected frequency bands. For the original 
OMP, however, the reconstruction errors are spread out in any frequency bands. 

When Sj is reconstructed at the j*'^ SU by OMP/Mod-OMP, a local occupancy estimate Oj needs to 
be obtained. For the k^'^ frequency band, 6^ (k) = 1 if at least one non-zero spike is detected. Otherwise 
dj (k) = 0: 

max||sj > Tf and i E Uk (14) 

where Tf is the threshold determined by the noise level. 

The fusion center then obtains 6 from all 6^. According to feature 4, the occupied frequency bands 
are common among all SUs. Though each SU estimates 6j independently, the original Oj are the same. 
Therefore, for frequency band k, we can use the L out of J rule for decision making [[TSl : 

J2d,{k) > L (15) 

j^l o(fe)=0 

In this correspondence, we set L = , where [•] is the smallest integer not less than •. This means that 
o{k) = 1 iff more than half of the SUs make the decision that frequency band k is occupied. Such rule 
can effectively reduce the number of wrong decisions, because each SU makes decisions independently, 
and the probability for L SUs to make the same wrong decisions is low. 
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In summary, the OMP/Mod-OMP DeF gets the estimate 6 by the following steps: 

1) Each SU reconstructs the spectrum Sj from by OMP/Mod-OMP. 

2) Each SU uses Sj to estimate an occupancy vector dj and sends it to the fusion center. 

3) The fusion center uses 6^ to determine 6 according to the L out of J rule. 

C. The Mod-SOMP for Data Fusion 

The proposed Mod-SOMP is modified upon the original SOMP. Mod-SOMP starts with Y, 0, r], u^. 
and C„. Mod-SOMP gets the estimate S as the following: 

1) Initialize: Set the iteration counter t = 0. Initialize the residual Res^*^ = Y, 0*^*^ = 0, the index 
set AW = Ci. If = 0, jump to Step 4. 

2) Obtain Z = [zi, Z2, ■ ■ ■ zj] using Equation dU 



^^Jt) (t) 



resl res2 ... res} 



3) Update the new residual matrix Res^*-* 

4) Increment t. Find the index A^*) that satisfies the following equation: 



using Equation [T2l 



= arg .jiiax^^ |(resf (16) 

5) If AW G Ufc c C2, let AW = A(*-i) U Ufc. If AW G c C3 or C„ = 0, let A^ = A(*-i) U A^*). 

6) Set the matrix of chosen atoms 0*^*^ = 0a{*)- 

7) Obtain Z using Equation [TT] and update the new residual Res*-*^ using Equation [121 

8) If t < m or Res^*-* ^ > V' return to Step 4. Otherwise, Z is the estimated signal matrix S with 
non-zero rows listed in A*^*\ 

Mod-SOMP becomes the original SOMP when C„ = and = 0. SOMP finds one non-zero index 
i in each iteration and treats all i independently. Mod-SOMP treats i with prior knowledge in the same 
way as Mod-OMP. Like Mod-OMP, Mod-SOMP has the same advantages over the original SOMP. 

When S is obtained, the fusion center first sums the columns of S to get s [[TSl : 

s(o = x:||s(*,j)|| (17) 

For the A;*'* frequency band, 6 (A;) = 1 if the following statement is true. Otherwise 6 (k) = 0. 

max ||s (z) II > Ty and i E Uk (18) 
In summary, the SOMP/Mod-SOMP DaF gets the estimate 6 by the following steps: 
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1) Each SU collects its own measurement vector and sends it to the fusion center, forming Y. 

2) The fusion center uses Y to reconstruct an estimate of the signal matrix S by SOMP/Mod-SOMP. 

3) The fusion center uses S to determine 6. 

IV. Simulation Results 

In this section, we will use the Monte Carlo simulation to compare the performances of the OMP DaF, 
Mod-OMP DaF, SOMP DeF, and the Mod-SOMP DeF. 500 simulations are conducted for each setup. In 
the simulations, there are = 41 frequency bands allocated to the PUs over a discrete spectrum with 
= 200. The frequency allocation information is defined by and C„. There are 2 in Ci, and 
||Ci||q = 5; 33 in C2 and ||C2||q = 165; 6 in C3 and HCsHq = 30. These numbers can be set 
according to the specific frequency allocation information in a specific spectrum range. 

In the simulation, a total of 10 are occupied randomly, with 2 in Ci, 2 in C2 and 6 in C3. The 
occupancy of the whole spectrum is around 10% (20 non-zero indices). Independent Rayleigh fading 
effects are applied at each SU, as well as Gaussian noise. 

We first simulate the J = 1, the single user case, with M = 70. In this setup, only OMP and Mod-OMP 
are compared. The R score vs. SNR plots are depicted in Fig. |2] From the R score we can see that the 
M = 70 is insufficient for OMP to get perfect correct decisions even when SNR is infinity. The R score 
for Mod-OMP, however, is increasing with SNR, and reaches 1 when SNR is infinity. 

Then, we simulate with J = 7 for all 4 CWSS algorithms with M = 70. The R score vs. SNR plots 
are depicted in Fig. |3] and Fig. |4l with M = 70 and M = 100, respectively. From Fig. |3l it can be seen 
that, the OMP DeF still cannot give perfect correct decisions when SNR is infinity. This is because none 
of the SU is able to provide correct decisions. The R scores of the other three algorithms, however, are 
able to reach 1 when SNR is infinity. The Mod-OMP DeF is better than the Mod-SOMP DaF, while the 
Mod-SOMP DaF is better than the SOMP DaF. This shows the advantage of using prior knowledge. 

From Fig. |4] We can see that with increased number of samples M, all algorithms have relatively better 
performance. The R scores for all 4 algorithms can reach 1 when SNR is infinity. The Mod-OMP DeF 
always has the highest R score. The Mod-SOMP DaF has higher R score than the OMP DeF when SNR 
is low, while lower R score than the OMP DeF when SNR is high. The SOMP DaF has lowest R score 
for all SNR. Generally speaking, the decision fusion algorithms (the OMP/Mod-OMP DeF) have higher 
R score than their data fusion counterparts (the SOMP/Mod-SOMP DaF). This is because decision fusion 
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makes independent decisions in each SU and has the L out of J rule to suppress wrong decisions. 

Notice that the SOMP DaF has very poor performance. We analyze one example of the decision statistics 
in Table HI] for M = 100, J = 7 and SNR = 15 dB. If all decisions are correct, a = 10, 6 = 0, c = 
and (i = 31. It can be seen that thanks to the L out of J rule, the decision fusion based algorithms 
make fewer wrong decisions on both occupied frequency bands and unoccupied frequency bands. The 
data fusion algorithms, however, cannot suppress wrong decisions on the unoccupied frequency bands. 
Then, it can be seen that both Mod-OMP DeF and Mod-SOMP DaF with the prior knowledge have less 
wrong decisions on the unoccupied frequency bands than their counterparts without the prior knowledge. 
This is because when the prior knowledge is applied, the algorithm will select one frequency band in each 
iteration, and the errors are restricted within the selected bands only. The algorithms without the prior 
knowledge, however, do not have such property and errors can be spread out in all frequency bands. 

A. Discussion 

A centralized CRN can choose between decision fusion and data fusion, depending on the required 
computation burden and communication burden. Decision fusion has high computation burden on each SU 
but low computation on the fusion center. Meanwhile, it has few communication burden because each SU 
only sends a binary vector o^ . Data fusion, however, has almost no computation burden on each SU but 
high computation burden on the fusion center. The communication burden is higher than that of decision 
fusion, because each SU needs to send the quantized sample vector y^. 

V. Conclusions 

In this correspondence, we have summarized four specific features of the CWSS and shown how to 
integrate those features in the algorithm design. With feature 1, we are able to use compressed sensing. 
With feature 2 and feature 3 as the prior knowledge, we are able to modify OMP and SOMP and 
proposed the Mod-OMP and Mod-SOMP. Then, with feature 4, we are able to integrate those algorithms 
in a centralized CRN. Four CWSS algorithms have been introduced, the OMP DeF, Mod-OMP DeF, 
SOMP DaF the Mod-SOMP DeF. Their decision performances in the CWSS have been compared by 
Monte Carlo simulations using the R score metric. From the simulation results, we have verified that 
the minimum number of measurements to achieve error free decisions can be reduced when the prior 
knowledge is applied. Then, we have verified that the decision fusion based OMP DeF and Mod-OMP 
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DeF have outstanding wrong decision suppression on both occupied frequency bands and unoccupied 
frequency bands. Finally, the Mod-OMP DeF and the Mod-SOMP DaF with the prior knowledge can 
restrict wrong decisions on unoccupied frequency bands. 
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TABLE I 

Properties of the Frequency Band Occupancy 





Time 


Frequency 


Category 1 


Always Occupied 


Fully Occupied 


Category 2 


Sometimes Occupied 


Fully Occupied 


Category 3 


Sometimes Occupied 


Partially Occupied 



Primary User 1 /' > j 



('t>) 



Secordary User 1 



"Primary User 3 



(CT>) 



I I ' 
I I I 



Secondary User 2 \^ Primary User 2 




Fusion Center 



Secondary User 3 



Serondary User 1 
Secontiary User 3 



M 



spectrum Observed at Different Secondary Users 



TABLE II 

Decision Statistics for All Algorithms 





Mod-OMP DaF 


OMP DaF 


Mod-SOMP DeF 


SOMP DeF 


a 


9.966 


9.982 


10 


10 


b 


0.034 


0.018 








c 


0.842 


2.906 


5.194 


11.944 


d 


30.158 


28.094 


25.806 


19.056 




Fig. 1. Illustration of a typical CRN coexisting with several PUs. The 
spectrum measured at all SUs have common non-zero entries. 



IVI = 70, J = 7 




Fig. 2. R Score vs. SNR for OMP and Mod-OMP in the single user 
WSS, M = 70. 



M = 100, J = 7 



Infinity 




6 9 12 

SNR (dB) 



infinity 



Fig. 3. R Score vs. SNR for the OMP DeF, Mod-OMP DeF, SOMP 
DaF and the Mod-SOMP DaF in the CWSS. M = 70, J = 7. 



Fig. 4. R Score vs. SNR for the OMP DeF, Mod-OMP DeF, SOMP 
DaF and the Mod-SOMP DaF in the CWSS. M = 100, J = 7. 



